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The choice of generator in the Similarity Renormalization Group (SRG) flow equation determines 
the evolution pattern of the Hamiltonian. The kinetic energy has been used in the generator for 
most prior applications to nuclear interactions, and other options have been largely unexplored. Here 
we show how variations of this standard choice can allow the evolution to proceed more efficiently 
without losing its advantages. 
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I. INTRODUCTION 

The Similarity Renormalization Group (SRG) uses a 
continuous series of unitary transformations to decouple 
high-momentum and low-momentum physics in an input 
Hamiltonian [TJ [2] . This decoupling means that expan- 
sions of physical observables generally become more con- 
vergent. The SRG can be implemented through a flow 
equation for the evolving Hamiltonian H S1 



dH s 
ds 



[r ls ,H s ] = [[G s ,H s \,H s ] 



(1) 



where s is a flow parameter |3] and the generator r] s 
is specified by the operator G s . With G s chosen to be 
the relative kinetic energy T re i, the SRG has been ap- 
plied successfully over the past few years to calculate nu- 
clear structure and reactions [4HT2] . However, different 
choices for G s will give rise to different patterns of evo- 
lution, which may be advantageous. In this paper, two 
alternatives to T re \ are evaluated for their effectiveness in 
decoupling and, in particular, for improvements in com- 
puting speed (see Fig. [I] for a representative example). 
Our tests are for realistic nucleon-nucleon (NN) inter- 
actions in two-body systems and for a one-dimensional 
model Hamiltonian applied to few-body bound states. 

We focus on novel generators that have G s as functions 
of T Te \. (Note: we can just as well consider the full kinetic 
energy T in our discussion, because the center-of-mass 
part commutes with the running Hamiltonian H s , so we 
will use T for convenience.) In particular, we explore the 
"inverse" G s operator 



1 + T/a 2 



g 1 : 



and the "exponential" G s given by 

G s = -a 2 e~ T ^ 2 = GT p • 



(2) 



(3) 



* |li.l287@osu.edu| 

t anderson@physics.ohio-state.edu 

* furnstahl.l@osu.edu 



1000 



100 



% 10 



0.1 



x— x G = T 



+ ■ - + G 
— ^ G 



exp 



J_ 



2.5 



eq 



3 

[fm 



3.5 



FIG. 1. (Color online) Computational time to evolve the Ar- 
gonne v±s 1 So potential 13 as a function of the final value 
of A eq (see text) for several generators, with a = 2 fm -1 for 
G\ nv and GT P - 



Each has a Taylor series that reduces to T (up to a con- 
stant, which drops out from the commutator) at low mo- 
mentum or when a is large. As such, the independent 
parameter a controls the separation of a low-energy re- 
gion where G s behaves as T and the potential is driven 
toward the diagonal, and a high-energy region where evo- 
lution is suppressed. This suppression can result in a 
significant computational speedup of the SRG evolution 
when compared to calculations with G s = T, while not 
impacting the advantageous properties of the evolution 
for low- momentum applications. 

The generators' suppression of running in unneeded 
parts of the Hamiltonian could mitigate the difficulties 
of evolving very large matrices for calculations of light 
atomic nuclei [11], opening the door to more tailored 
SRG generators and more effective evolution of three- 
and eventually four-body interactions. Even at the two- 
body level there are problems when trying to evolve to 
large values of the flow parameter s. A recent example is 
a study of SRG decoupling with large-cutoff effective field 



2 



theory (EFT) potentials, which require evolution beyond 
the range usually considered [14]. In doing so, the SRG 
differential equations can become extremely stiff and take 
a prohibitively long time (weeks on a single processor) to 
evolve. This problem has hindered exploratory studies 
into issues such as what happens when a chiral EFT is 
evolved to the regime of pionless EFT. 

In Sec. [TTJ we give representative results for applica- 
tions to two-body systems, including an analysis of the 
flow pattern. These results are extended to few-body sys- 
tems in Sec. |III| using a model one-dimensional Hamilto- 
nian that has proved useful in past applications p~5j [16] . 
We summarize and outline future studies in Sec. IIVI 



II. TWO-NUCLEON SYSTEMS 

In this section we give representative results for evolv- 
ing realistic nucleon-nucleon potentials using the novel 
generators from Eqs. Q and (|3| in comparison to the 
usual choice of G s = T. 



A. Performance 

The key advantage of the generators that we high- 
light is the improvement in computational performance. 
For example, the time needed to evolve the Argonne v\% 
1 So potential [13] to equivalent levels of decoupling with 
several generators is plotted in Fig. [I] The parameter 
A = 1/s 1 / 4 , which has dimensions of a momentum, has 
been used to identify the momentum decoupling scale. 
However, different generators will evolve a given poten- 
tial at different rates, so comparing results with the same 
definition of A can be misleading. Therefore, we iden- 
tify an "equivalent" A eq for each generator that equal- 
izes the degree of decoupling compared to G s = T (for 
which A eq = A by definition); the details are described in 
Sec. iHBl 

The value of a will also have an impact, as dis- 
cussed below; in Fig. [I] we use the intermediate value 
<j = 2fm _1 . With this choice, there is nearly an or- 
der of magnitude difference in the time to evolve the 
Argonne vi$ potential with GJ xp and G x s nv compared to 



4fm 1 and two orders of magnitude 
Note that nuclear interactions typi- 



G s — T at A eq 
by A eq = 2fm" 
cally have been evolved for nuclear structure studies in 
the range A = 1.5-2.2 fm _1 . The speed gains will de- 
pend on the initial potential and can be much less for the 
evolution of softer initial potentials; e.g., for the N 3 LO 
500 MeV chiral EFT potential of Ref. [17 , evolving with 
GJ xp to A eq = 2fm _1 is about 1.5 times as fast as with 
G s =T and about 3 times as fast to A eq = 1.5 fm - . 

The numerical solution of the SRG evolution equa- 
tions requires repeated dense matrix-matrix multiplica- 
tions. As a consequence, the evolution has been car- 
ried out on shared- memory computer architectures. Re- 
cent calculations using the SRG with many-body forces 



are approaching the limits of what is practical to evolve 
in memory on a single node because of the size of the 
model space needed [11]. A distributed scheme to solve 
the equations would permit larger model spaces to be 
utilized; however, the dense matrix multiplication would 
then be limited by internode communication times. The 
reduced number of operations required by novel genera- 
tors might help make such a scheme possible. 



B. Decoupling and A eq 

To validate the apparent computational advantages of 
these generators, one must confirm that the decoupling 
characteristics of the G s = T generator are also repro- 
duced, so that calculations of physical observables also 
become more convergent. However, if we evolve to the 
same A, the degree of decoupling for identical initial po- 
tentials differs for GJ xp and G x s nv compared to G s = T. 
These differences are evident in the deviations of 1 So 
phase shifts calculated from the evolved potentials using 
G® xp and G s = T from the unevolved potential shown in 
Fig. [2] for several different A values (only G® xp is shown; 
G™ behaves similarly). If the full potentials were used, 
the phase shifts would agree - up to numerical preci- 
sion - with those from the initial potential, because the 
evolution in all cases is unitary. However, the degree of 
decoupling for a given value of A can be made manifest 
by first cutting off the potential (that is, setting its ma- 
trix elements to zero) above some value of k and then 
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FIG. 2. (Color online) Deviation of phase shifts calculated 
with the Argonne vis 1 So potential [13] evolved with G s —T 
and G| xp (a = 2 fm -1 ) to various A values and then trun- 
cated at fccut = 2 fm -1 to test decoupling. Phase shifts from 
untruncated potentials agree precisely with those from the 
initial potential. 
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FIG. 3. (Color online) Contour plots showing the evolution via Eq. |l]) at several values of A eq starting from the momentum- 
space Argonne vis 1 So potential 13 using (a) G s = T, (b) GT from Eq. (2) with a = 2 fm -1 , and (c) Gf p from Eq. Q with 
cr = 2fm" 1 . 



calculating the phase shifts. In Fig. [2j for illustration, 
we choose the cutoff value fc cut to be 2fm _1 . The signa- 
ture of decoupling is that the phase shifts agree at lower 
energies and only deviate close to and above the cutoff. 
This is typically observed when the potential is evolved 
so that A is less than the cut momentum p~8j [19] . 

This behavior provides us with a pragmatic way to 
define A eq : identifying it with the decoupling behavior 
that is found for G s = T for a given A. We use the results 
in Fig. [2] to illustrate the procedure. The continuous 
curve shows the deviation of phase shifts for the potential 
evolved with G s = T to A = 2 fm -1 and cut at k = 
2 fm -1 . As expected from decoupling, the deviation is 
small up to roughly fc cut . We use this level of agreement 
as the criterion for identifying equivalent A's for other 
generators. That is, a potential is evolved to a series of 
A values with a different G s and then cut at fc cut after 
evolution. The phase shifts are then compared to the 
level of decoupling observed for G s = T at k cut = A. The 
approximate point in the novel G s evolution for which the 



cut phase shifts agree is equated with A eq . Consider the 
phase shifts for the potential evolved with GJ xp to several 
values of A and cut at k = 2 fm -1 , as shown in Fig. [2J 
For G® xp , A evolved to 1.5 fm -1 gives a similar degree of 
decoupling as A evolved to 2 fm -1 with G s = T. As a 



result, we define A = 1.5 fm 



to be the A eq = 2 fm 



for 



G e s xp . We use A eq for most comparisons in the following 
discussion. 



C. Flow Analysis 

Having chosen a working definition for the decoupling 
scale of evolution with our novel generators, we take a 
closer look at the properties of the evolved potentials. In 
particular, we would like to see how a potential flows with 
evolution using these generators compared to G s = T and 
to understand how the choice of a affects this flow. 

In Fig. [3] we compare the evolution pattern of the two- 
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FIG. 4. (Color online) Diagonal (a) and off-diagonal (b) momentum-space matrix elements for evolved Argonne vis 1 So 
potential [13] with G s = T and generators Gf p and G!, nv at A = 2fm _1 with a = 2 fin -1 . 
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FIG. 5. (Color online) Diagonal (a) and off-diagonal (b) momentum-space matrix elements for evolved Argonne v\% 1 So 
potential [13] with G s = T and G| xp with different values of a each evolved to A = 2fm _1 . The lines with a = 3 fm -1 are 
indistinguishable from the ones with G s = T. 



body potential in the x So channel with different gen- 
erators. Each frame is a representation of the poten- 
tial matrices in momentum space; where the matrix is 
zero, there is no coupling between momentum compo- 
nents. The initial potential in all cases is Argonne vig [T3] 
and the value of a is taken to be 2 fm -1 . Note that at 
A = 4fm _1 , the first matrix plotted here, there is already 
significant evolution. As the potential is evolved, its high 
and low momentum components become increasingly de- 
coupled, as expected. It is evident that the evolved po- 
tentials are similar (but not identical) in the region where 
k 2 ,k f2 < a 2 . When A < a we find that a roughly defines 



the low-momentum region where these generators behave 
as G s = T. The minor differences in this low-momentum 
region can be attributed to the fact that the point in 
evolution, A, needed for these generators to reach the 
corresponding A eq occurs when A < A eq . At higher mo- 
menta, novel generator evolution is suppressed relative 
to G s = T. The patterns here are characteristic of the 
particular generator and are similar for other potentials 
and in other channels. 

This nature of the evolution is further illustrated by 
the plots in Figs. [4] and |5j which show a detailed view 
of the diagonal and off-diagonal values of the matrices 



5 



k' 2 (fm~ 2 ) 





^ 3 
E 

H 

^ 6 
9 



J 



(a) 



1 0.5 

0(fm) 
-0.5 



k' 2 (fm 2 ) 



k' 2 (fm 2 ) 



k' 2 (fm 2 ) 



^ R ( 
\j \j \j \. 


) ^ 




) 3 6 


9 










G ="M 

S 


G inv 1 

S 




G exp l 

S 


_ 



(b) 



1 0.5 



0(fm) 



-0.5 



FIG. 6. (Color online) (a) Unevolved N 3 LO 500 MeV potential . (b) Contour plot showing the evolved N 3 LO 1 S potential 
with different generators to A eq = 2fm _1 . For GT P and , a = 2fm _1 . 
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FIG. 7. (Color online) Color contour plots of the first term 
(excluding the factor due to V s ) on the right side of Eq. ^ 
for G s = T (left), Gf" (middle), and G? p (right). The last 
two use a — 2 fm -1 . 
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FIG. 8. (Color online) Color contour plots of the first terms 
(excluding the factor due to V s ) on the right side of Eq. ^ 
for GT P with a = 1,2, and 3fm _1 from left to right. 



at low momentum for each generator applied to the Ar- 
gonne v±g 1 $o potential . The flow parameter was run to 
A = 2fm _1 with the generators here to distinguish be- 
tween ambiguities caused by using A eq . In Fig. HI the val- 
ues from different generators agree quite well for k < a, 
where a ~ A. Moreover, we see that the minor "pincush- 
ion" effect seen in Fig. [3] for these generators relative to 
G s = T is indeed an artifact of using A eq , as each curve 
falls to zero at approximately the same time. In Fig. [5] 
we focus on GJ xp versus G s — T and the fact that they 
differ for a < A. However, as a — >• A it is evident that 
the evolution of the novel generator becomes increasingly 
similar to G s = T in the low- momentum region k < A. 
For larger values of a, these generators become indistin- 
guishable from G s = T. 



As a demonstration of similar behavior for different NN 
potentials, an evolved potential with different generators 
in the % channel of the N 3 LO 500 MeV potential [17] 
is shown in Fig. [6j Note that the initial potential has 
significantly less coupling at high momentum compared 
to Argonne v\g [8 . As a result, there is correspond- 
ingly less improvement in evolution speed. However, the 
general features of the evolution patterns with different 
generators seen with Argonne vig are also seen for the 
N 3 LO potential. 

To better understand the evolution process, we need 
to look further into the flow equation itself. Evaluat- 
ing Eq. in a two-body partial-wave momentum space 
basis with G s = T yields 

dV s (k,k>) = _ {k 2_ k ,2 ) 2y s{Kk , ) 

+ - / q 2 dq(k 2 + k' 2 -2q 2 )V s (k,q)V s (q,k f ) .(4) 
* Jo 

In the far off-diagonal region, the first term dominates 
(this is true for the ordinary range of A but is modified 
when A is comparable with the binding momentum of a 
bound state). This implies that each off-diagonal matrix 
element is driven to zero as 

V s {kX) V s=0 (k,k')e-^ k2 - k ' 2 ^ . (5) 
For G s = f(T), these results are modified to 

dVs ^ k>) = -(fc 2 - fc' 2 ) (/(fc 2 ) - f(k"))V s (k, k') 

+- r^(/(fc 2 )+/(fc' 2 )-2/( 9 2 )) 

K Jo 

xV s (k,q)V s (q,k f ) (6) 

and 

Vs (k, fc') V s=0 (k, fc') e -^ 2 -fc' 2 )(/(* 2 )-/(fc' 2 )) . ( 7 ) 

The difference in the exponents of Eqs. ([5| and Q for 
k ~ A leads to A eq < A. 

The first term of the flow equation (excluding the fac- 
tor due to the potential) for each of our generators is 
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FIG. 9. (Color online) Contour plot showing the evolution to A eq = 2fm 1 via Eq. |T]) starting from the momentum- space 
Argonne vis 1 So potential [13] using G| xp from Eq. ^ with different values of a. 
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FIG. 10. (Color online) Computational time to evolve the 
Argonne v±s 1 So potential [13] with G s — T and G| xp to 
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shown as a contour plot in Figs. [7] and [8] In Fig. [7] we 
see that T works uniformly (in k*) on the entire region 
of the potential. With G 1 /^, there is much less evolution 
close to the diagonal. The plot for G® xp is similar, but 
exhibits even less evolution in the middle region. 

It is evident here that the novel generators will re- 
sult in less evolution at high momenta, as we have seen, 
and that this should be a generic result for other partial 
waves and for higher-body evolution, because it depends 
on only kinetic energy differences. We also see how the 
value of a controls the degree to which the operator G s is 
similar to T. This is illustrated in Fig. [SJ If a = 1 fm -1 , 
only the edges of the potential are modified; the shape 
is completely different from T. For a = 3 fm -1 , there is 
the thinnest band on the diagonal, which is closest to T. 
At very large <r, there is a transition to T. In the plots 
of Fig. H we see how the final evolved flow is affected by 
differing choices of a. 



The limited evolution at high momenta seen in these 
generators suggests that the time to evolve to a given 
decoupling parameter A eq should be less for G x s nv or G^ xp 
than for G s = T. The dramatic drop in evolution time 
seen in Fig. [I] for the generators makes it apparent that 
they are more efficient. However, we can also look at how 
the choice of a affects their performance. This is shown 



in Fig. 10 for G| xp , where we see that the time spent 
evolving to A eq = 2 fm -1 decreases as a decreases. As a 



becomes smaller, the evolution at high momentum is in- 
creasingly limited, which we correlate with improvement 
in computation time; as a increases, the evolution time 
approaches that of T (as does the flow). Note that for 
the Argonne vis potential a large a is needed before G s 
is effectively equal to T. 

In practical applications, one might optimize the trade- 
off between decoupling and computational speedup by 
choosing a to be approximately equal to, or slightly 
greater than the A corresponding to A eq . Then the de- 
coupling properties of G s = T are preserved in the low 
momentum region of interest while still enhancing the 
computational performance of the evolution by limiting 
the evolution at high momentum. However, this prescrip- 
tion has not yet been tested in detail. 



III. FEW-BODY TESTS IN A 
ONE-DIMENSIONAL MODEL 



The effects on matrices in a momentum basis demon- 
strated in the last section are generic and so should carry 
over to alternative bases and to higher-body forces. Cal- 
culations for realistic three-dimensional few-body sys- 
tems are not yet available, but we can test the generators 
in a one-dimensional model of bosons that has proven to 
accurately predict the evolution of three-dimensional few- 
body forces [15]. The model we use was originally intro- 
duced in Ref. [20 as a sum of two Gaussians to simulate 
repulsive short-range and attractive mid-range nucleon- 
nucleon two-body potentials. It is written in coordinate 
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Basis 


Momentum 


Oscillator 


System 


A = 2 


A = 2 


A = 3 


A = 4 


a = 2 


3.3 


3.3 


3.8 


4.1 


a = 3 


2.6 


2.6 


2.7 


2.8 



TABLE I. Speed up in model A-particle one-dimensional 
(ID) oscillator basis for evolution to X eq ~ 3, comparing the 



ratio of the time to evolve G s 
for are very similar. 



T versus Gf p . The results 



space as 



^( 2 )(x) 



v 2 



and in momentum space as 



(8) 



(p-p')M/8 



V 2 



2^^ 



2 e -(p-p X ) 2 ^/8 



(9) 

where Vi 12, = 0.2, F 2 = -12, and a 2 = 0.8 (see 
Ref. [20 for discussion of units). Also, in calculations 
with an initial three-body potential, a regulated contact 
interaction is used. This is written as 

V^(p,q,p',q') = c E f A (p,q)f A ( P ',q'), (10) 

where ce corresponds to the strength of the interaction 
and 



h{p,q) = e-^ + ^l^\ 



(11) 



The regulator cutoff A = 2, and the sharpness of the fall 
off is set to n = 4. We have chosen these parameters for 
comparison with Ref. [T5] . 



A. Performance 

Our first test is to confirm that the enhanced com- 
putational performance characteristics of these genera- 
tors are maintained in the few-body basis. Table [I] shows 
the speedup obtained using the GJ xp generator with the 
model potential described previously. The performance 
of the G x s nv generator is very similar to the G^ xp results 
quoted here. Results are reported at an optimal A eq ; we 
do not use a range of A values here because of the par- 
ticularities of the oscillator basis convergence properties, 
which is discussed in more detail in what follows. 

We find that the evolution of the two-body force in 
two-, three-, and four-particle systems with novel gener- 
ators are all 2.5-4 times faster than the evolution with 
G s = T (to the same degree of decoupling). The per- 
formance enhancement is relatively basis independent, 
with the speedup for the momentum and oscillator bases 
roughly equivalent. While the speedup improvement is 
much smaller than that found for Argonne t^g, we do not 
expect a direct correspondence. The important point is 



that the speedup in the A=2 particle system serves as a 
good predictor of the speedup in the A=3 and A=4 sys- 
tems. Thus, one might expect a similar improvement in 
few-body oscillator basis calculations with Argonne v\% 
as found previously for Argonne vig in the two-particle 
partial-wave momentum basis. This will be significant 
as we move to novel generator calculations in realistic 
three-dimensional systems. 



B. Decoupling 

The measure of performance using these generators de- 
pends explicitly on their decoupling properties in the few- 
body harmonic oscillator basis relative to G s = T. Ul- 
timately, we find the level of decoupling obtained with 
G s = T to be matched by these generators. 

However, the convergence properties of the oscillator 
basis with respect to SRG evolution, and consequently 
the issue of selecting an appropriate A eq in this basis, is 
more complicated than for the momentum basis. The 
convergence of observables depends on a balance of the 
ultraviolet (UV) and infrared (IR) cutoffs intrinsic to the 
choice of a particular oscillator basis. These cutoffs are 
given by [TT] 



A uv ~ y/mN max Ml 



and 



A 



IR 



rnhVt 



(12) 



(13) 



where ft is the oscillator frequency and Af max is the max- 
imum number of total oscillator excitations in the ba- 
sis. Thus, a cutoff in oscillator basis states results in 
two approximate cutoffs in momentum space. However, 
the SRG, using the generators considered here, provides 
only a means to effectively lower the UV cutoff (by de- 
coupling high- and low-momentum degrees of freedom in 
the Hamiltonian). As such, convergence is not monoton- 
ically improved with respect to evolution in A. 

As a measure of the decoupling, we plot the binding 
energy of the lowest energy state for an A-particle system 
with respect to iV cu t for evolutions of the initial Hamil- 
tonian to various A (this procedure was carried out for 
G s = T using this model in [E]). The actual calculation 
is carried out by evolving the model interaction to A in 
an initial basis large enough so that the binding energy 
is well converged. The Hamiltonian is then truncated at 
N cut and the binding energies are calculated in the re- 
duced basis. The value of N cut refers to the number of 
oscillator excitations in the basis and is the oscillator ba- 
sis equivalent of the fc cut parameter in momentum space, 
as used in Sec. [Hj Again, each N cut here corresponds to 
a rough Auv and Air truncation in momentum space. 

Results are shown in Fig. 11 for the A = 2, A= 3, and 
A=4 particle systems using GJ xp with <j = 3 for selected 
values of A. The signature of decoupling is the improved 
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FIG. 11. (Color online) Decoupling using G e * p with a = 3 for (a) A = 2, (b) A = 3, and (c) A = 4. The initial two-body-only 
potential is evolved to each A shown in a basis with A^ max = 40. Matrix elements of the potential are set to zero if one or both 
states have N > N cut and the resulting Hamiltonian is diagonalized to obtain the ground-state energies plotted. 




FIG. 12. (Color online) The lowest bound-state energy for a four-particle system as a function of A evolved using the SRG 
Gl xp generator and an A^ max = 40 basis with an initial three-body potential with strength (a) ce = 0.0, (b) ce = —0.05, and (c) 
ce = 0.05. The (blue) curves with diamonds include the full evolution of the Hamiltonian while the (black) curves with circles 
use the two-body potential evolved in the two-particle system and the (red) curves with squares use the two- and three-body 
potentials evolved in the three-particle system. Evolution with is almost indistinguishable. 



convergence of the binding energy at smaller N cut with 
respect to SRG evolution. As the interaction is evolved, 
the degree of decoupling gets better. This is true only up 
to some value of A, however, at which point the degree of 
decoupling starts to get worse. It is the latter behavior, 
introduced by the IR cutoff of the oscillator basis, that 
complicates our efforts to choose a A eq because a one-to- 
one correspondence with A is no longer clear. 

Nevertheless, a practical choice can be made by equat- 
ing A with A eq when decoupling is found to be optimal 
in the G s = T and novel generator evolutions. The op- 
timal levels of evolution happen to coincide at A « 3 for 
G s = T and these generators. Thus, the speedup results 
in Table [I] were reported at A eq = 3. Moreover, given 
these values for A and A eq , we have chosen a = 3 for 



most of the model space calculations in this section. 

One may note that differences do exist between the 
G s = T and novel generator decoupling results, particu- 
larly at low N cut [15]. However, the level at which any of 
the generators become well converged to the exact results 
with respect to N cut are effectively the same. Thus, it is 
reasonable to make the comparisons we have done here 
to determine A eQ . 



C. Induced Many-Body Forces 

In general, the evolution of an interaction via the SRG 
leads to induced many-body forces. This is evident if we 
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examine the second quantized form of the Hamiltonian 



H 



Tija\cij 



Vijkia\a^aia k 



(14) 



where the dots indicate that higher-body forces may also 
be present in the initial interaction. When the SRG 
commutators in Eq. are performed, one can see that 
many-body forces will be induced. These induced forces 
could pose a serious problem for few-body calculations 
because at some point we must truncate the model space 
in numerical calculations, which will alter the predicted 
value of observables. These can be controlled, however, 
if there is a hierarchy of many-body forces so that suc- 
cessively larger many-body components are suppressed, 
and one can include at most one or two induced pieces to 
obtain well converged results. This has been found to be 
the case for G s = T and needs to hold for any practical 
alternative generators. 

A measure of the induced many-body forces (which has 
been used in previous studies [6] [TTJ [15]) is to calculate 
the ground-state energies of the A = 3 and A=4 parti- 
cle systems starting with a two-body interaction and to 
examine how this energy changes with and without the 
induced three- and four-body components, as a function 
of the evolution parameter A. We do this in Fig. 12 with 
a plot of the ground-state energy for the four-particle sys- 
tem with the initial two-body interaction embedded and 
evolved in the A = 2, A=3, and A=4 body bases with 
the G| xp generator. The results for G x s nv are virtually 
indistinguishable. The curves show that the hierarchy of 
induced many-body forces is preserved for these genera- 
tors just as with G s = T [15 (see, however, Ref. [2T]). 
This hierarchy also holds for calculations with an initial 
three-body force and in the A = 3 particle system. 

In summary, the model results suggest that the ad- 
vantageous features of SRG evolution with G s = T 
can be maintained with the added computational per- 
formance of these generators when applied to realistic 
three-dimensional few-body calculations. 



IV. SUMMARY 

In this work, novel generators for the SRG that are 
functions of the kinetic energy operator T with an ad- 
justable scale parameter a were tested. We found that 
functions that reduce to T for basis states with kinetic 
energy less than a preserved the good features of T, such 
as decoupling, but efficiently suppressed evolution for 
higher kinetic energies and thereby took much less time 
to evolve. Specific examples were considered, but other 
choices with a Taylor expansion starting with T should 
give comparable results. Their action was understood us- 
ing a simple analysis of how the generators directly affect 
regions of high and low momentum. If a is large enough, 



the generators become equivalent to T. It is important 
to note that not only the two-body properties of T were 
preserved by these generators, but so were its characteris- 
tics in a few-body model space. This includes decoupling 
and the hierarchy of induced many-body forces, which is 
critical for applications to larger systems of particles. 

The generators allow us to evolve potentials to much 
smaller values of A than previously feasible. This should 
enable us to explore the transition between pionful and 
pionless regions of EFT potentials and further test the 
observations of Glazek and Perry about evolving past a 
bound state [22] . The original choice for G s advocated by 
Wegner and collaborators [2] [3] and applied extensively 
in condensed matter is the diagonal component of the 
interaction, G s = Hdi& g (s), 



(i\H diag (s)\j) 



(i\H(s)\j) ifi=j 







otherwise. 



(15) 



In Ref. [22], it was observed that when evolving a sim- 
ple model past a bound state the Wegner evolution with 
#diag will decouple the bound state by leaving it as a 
S function on the diagonal of the Hamiltonian. In con- 
trast, with G s = T the bound states remained coupled 
to low momentum and were pushed to the lowest mo- 
mentum part of the matrix. This behavior was explored 
in Ref. [14] for leading-order, large-cutoff EFT poten- 
tials featuring deeply bound spurious states. However, 
it has not been studied for the physical deuteron state, 
which requires evolving well below A = 1 fm _1 . This is 
now easily possible with the replacement of i^diag f° r T 
in Eqs. Q and (J3|, although there are as- yet-unsolved 
complications from the discretization of the momentum 
basis. 

The most important next step for the novel generators 
is to apply them to evolve realistic few-body potentials, 
where speeding up the evolution is desirable because of 
the large sizes of the matrices involved. The generators 
can be applied directly to few-particle bases using the 
method described in Refs. [6j [11] ITS] , Calculations in 
a one-dimensional model performed here imply that the 
speed up carries over to three-body forces and could have 
a significant impact in making realistic calculations with 
additional induced many-body forces feasible. 
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